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The operation of atomic clocks is generally carried out at room temperature, whereas the definition 
of the second refers to the clock transition in an atom at absolute zero. This implies that the clock 
transition frequency should be corrected in practice for the effect of finite temperature of which 
the leading contributor is the blackbody radiation (BBR) shift. In the present work, we used 
configuration interaction + coupled-cluster method to evaluate polarizabilities of the 6s 2 1 So and 
6s6p 3 P states of T1+ ion; we find q ( 1 5'o) = 19.6 a.u. and a ( 3 Po) = 21.4 a.m. The resulting BBR 
shift of the 6s6p 3 P - 6s 2 1 5 T1+ transition at 300 K is A^bbr = -0.0157(16) Hz. This result 
demonstrates that near cancelation of the 1 So and 3 Po state polarizabilities in divalent B + , Al + , 
In ions of group IIIB [Safronova et at, PRL 107, 143006 (2011)] continues for much heavier T1+, 
leading to anomalously small BBR shift for this system. This calculation demonstrates that the 
BBR contribution to the fractional frequency uncertainty of the Tl + frequency standard at 300 K 
is 1 x 10 -18 . We find that Tl + has the smallest fractional BBR shift among all present or proposed 
frequency standards with the exception of Al + . 

PACS numbers: 06.30.Ft, 31.15.ac, 31.15.ap, 31. 15. am 



I. INTRODUCTION 

Recent advances in atomic and optical physics have 
led to unprecedented improvements in the accuracy of 
optical frequency standards that are essential for many 
applications including measurements of the fundamental 
constants and search of their variation with time, testing 
of physics postulates, inertial navigation, magnetomctry, 
tracking of deep-space probes, and others [if. An optical 
clock with a record low fractional frequency uncertainty 
of 8.6 x 10 -18 , based on quantum logic spectroscopy of 
an Al + ion was demonstrated in 2010 [2j. 

Any definition of the second should be based on a clock 
decoupled from its particular environment. Thermal fluc- 
tuations of the electromagnetic field, i.e. blackbody ra- 
diation (BBR) , are pervasive and can only be suppressed 
by cooling the clock. The BBR at any non-zero tempera- 
ture induces small shifts in atomic energy levels through 
the AC Stark effect. The operation of atomic clocks is 
generally carried out at room temperature and the clock 
transition frequency should be corrected in practice for 
the BBR shift. Experimental measurements of the BBR 
shifts are sufficiently difficult that no direct measurement 
has yet been reported for optical frequency standards. At 
room temperature, the BBR shift of a clock transition 
turns out to make one of the largest irreducible contribu- 
tions to the uncertainty budget of optical atomic clocks 
Q. The present status of the theoretical and experi- 
mental determinations of the BBR shifts in all frequency 
standards was recently reviewed in @, |3| ■ 

The BBR frequency shift of a clock transition can be 
related to the difference of the static electric-dipole po- 
larizabilities between the two clock states Recent 
work [6] demonstrated that the polarizabilities of ground 



ns 2 1 Sq and metastable nsnp 3 P states are nearly equal 
to each other in B + , Al + , and In + , all of which are group 
IIIB ions. As a result, these three ions have anomalously 
small BBR shifts of the ns 2 1 Sq — nsnp 3 Po clock tran- 
sitions. The fractional BBR shifts for these ions are at 
least 10 times smaller than those of any other present or 
proposed optical frequency standards at the same tem- 
perature, and are less than 0.3% of the Sr clock shift. 

Optical frequency standard based on 204^1+ 
6s 2 1 So nip — - 6s6p 3 Po nip> = transition 
was proposed in Ref. Q. The radioactive isotope of 
204 T1 has a half-life of 3.78 years, a spin of 2, and a very 
small magnetic moment of 0.0908 nuclear magnetons 
making it ideal object for very high-resolution laser 
spectroscopy Q- Because of its small nuclear magnetic 
moment the natural linewidth of the clock transition in 
204rpi+ > g ex p ec t ec i to be orders of magnitude smaller 
than estimated for stable Tl isotopes The BBR 

in this frequency standard have not been previously 
estimated. Since three group IIIB ions exhibit very 
small BBR shifts, it is very interesting to evaluate if this 
trend holds for much heavier Tl + . 

The BBR frequency shift of the clock transition can 
be related to the difference of the static electric-dipole 
polarizabilities between the clock states, Aa , by @ 

A^bbr = -~(831.9 V/m) 2 (^-^ Aa (l + „), (1) 

where 77 is a small dynamic correction due to the fre- 
quency distribution and only the electric-dipole transi- 
tion part of the contribution is considered. The Ml and 
E2 contributions have been estimated for Al + and found 
to be negligible l 6j]. Therefore, the calculation of the BBR 
shift reduces to accurate calculation of the static polar- 
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izabilities of the clock states and dynamic correction rj. 

In this work, we evaluate polarizabilities of the 6s 2 
and 6s6p 3 Po states in Tl + , corresponding BBR shift and 
its uncertainty. Dynamic correction to the BBR shift is 
evaluated and found negligible. We also calculate a num- 
ber of electric-dipole matrix elements in Tl + for transi- 
tions between low-lying levels. We note that our cal- 
culation of all of these properties is independent on the 
particular isotope number well within the quoted level of 
precision. Therefore, all these results apply to any T1+ 
isotope. 

II. METHOD 

Correlation corrections between a few valence electrons 
can be accurately treated by the configuration interac- 
tion (CI) method. Since the valence-valence correla- 
tions are very large, the CI method provides better de- 
scription of these correlations than the perturbative ap- 
proaches. However, excitations of the core [Is 2 , 5d 10 ] 
electrons can not be directly included in the CI approach 
due to enormous size of such problem. An elegant ap- 
proach to the inclusion of the core-valence correlations 
within the CI framework was developed in Q , where core- 
valence correlations were incorporated into the CI by con- 
structing an effective Hamiltonian using the second-order 
many-body perturbation theory (CI+MBPT). Recently, 
we have developed the relativistic Cl+all-order method 
combining CI with coupled-cluster (CC) approach. 
This method, first suggested in [l(|, was successfully ap- 
plied to the calculation of divalent atom properties in 
Refs. @, Q. The coupled-cluster method used here is 
known to describe the core-core and core-valence corre- 
lations very well as demonstrated by its grea t success 
in predicting alkali-metal atom properties [Tl|. There- 
fore, combination of the CI and all-order coupled-cluster 
methods allows to account for all dominant correlations 
to all orders. To evaluate uncertainty of our results, we 
use all three of the approaches and compare the results of 
the CI, CI+MBPT, and Cl+all-order calculations. We 
refer the reader to Refs. [HH,® for the description of the 
methods and outline only main points of the calculations 
below. 

We start with solving Dirac-Fock (DF) equations 

H ip c = e c ip c , 

where Hq is the relativistic DF Hamiltonian [H, [9| and ip c 
and e c are single-electron wave functions and energies. 
The self-consistent calculations were performed for the 
[Is 2 , 5cZ 10 ] closed core and the 6s, 7s, 6p, 7p, and 6d 
orbitals. We have constructed the B-spline basis set con- 
sisting of N = 35 orbitals for each of the s, Pi/%, P3/2, ■■■ 
partial waves up to I < 5; core, 6s, 7s, 6p, 7p, and 6d 
orbitals were replaced by the exact DF functions for in- 
creased accuracy. The basis set is formed in a spher- 
ical cavity with radius 60 a.u. The CI space is effec- 
tively complete and includes 20sp and 21dfg orbitals. All 



MBPT and all-order terms were summed over the entire 
N = 35, I < 5 basis set. 

The multiparticle relativistic equation for three valence 
electrons is solved within the CI framework [l2[ to find 
the wave functions and the low-lying energy levels: 

H cS (E n )<S> n = E n $ n . 

The effective Hamiltonian is defined as 

H cS {E) =H FC + X(E), 

where H-pc 18 t ne Hamiltonian in the frozen-core approx- 
imation and the energy-dependent operator T,(E) takes 
into account virtual core excitations. The £(£7) part of 
the effective Hamiltonian is constructed using the second- 
order perturbation theory in the CI+MBPT approach 
Q and linearized coupled-cluster single-double method 
in the Cl+all-order approach [9|. The T,(E) = in the 
pure CI calculation. Construction of the effective Hamil- 
tonian in CI+MBPT and Cl+all-order approximations is 
described in detail in Refs. @,H|. 

III. RESULTS 

Comparison of the energy levels (in cm -1 ) obtained 
in the CI, CI+MBPT, and Cl+all-order approximations 
with experimental values [HHHI is given in Table Q] Cor- 
responding relative differences of these three calculations 
with experiment are given in the last three columns in 
%. Two-electron binding energies are given in the first 
row of Table HI energies in other rows are counted from 
the ground state. We also observed significant, by a fac- 
tor of 4 or better, improvement in the precision of the 
energy levels with Cl+all-order method in comparison 
with the CI+MBPT one. For example, CI+MBPT value 
for the two-electron binding energy differs from experi- 
ment by 1.8%, while our all-order value differs from the 
experiment by only 0.4% (see line one of Table H]). The 
experimental value of the two-electron binding energy 
is obtained as the sum of the Tl + and Tl 2+ ionization 
limitsgiven in [lj], 164765(5) cm" 1 and 240600 cm" 1 . 
Ref. [13[ notes that the ionization limit for Tl 2+ , derived 
from the first 3 members of the 2 S series was shifted by 
300 cm _1 to give effective quantum number for 5g 2 G 
that is nearly hydrogenic. Therefore, there is some uncer- 
tainty (<■ 0.1%) associated with the two-electron binding 
energy in Tl + . 

We also compared the transition energies between the 
6s6p 3 Po level and 4 levels relevant to the calculation of 
the 6s6p 3 -Po polarizability. These values, calculated in 
the Cl+all-order approximation are compared with ex- 
periment in Table [XT] We find that these transition en- 
ergies are substantially more accurate than the energy 
levels counted from the ground state listed in Table HI 

In the present calculation, the Tl + scalar polarizabil- 
ity cto is separated into a valence polarizability , ionic 
core polarizability a c , and a small term a vc ( that mod- 
ifies ionic core polarizability due to the presence of two 
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TABLE I: Comparison between experimental [LsllTil] and the- 
oretical energy levels in cm -1 . Two-electron binding energies 
are given in the first row, energies in other rows are counted 
from the ground state. Results of the CI, CI+MBPT, and 
Cl+all-order calculations are given in columns labeled CI, 
MBPT, and All. Corresponding relative differences of these 
three calculations with experiment are given in the last three 
columns in %. 



State 


Expt. 


CI 


MBPT 


All 


Differences (%) 
CI MBPT All 


6s' 1 


So 


405365 


376102 


412676 


407125 


-7 


1.8 


0.4 


6s7s 


3 Sl 


105229 


92945 


108031 


106028 


-12 


2.7 


0.8 


6s7s 


's 


108000 


96304 


110845 


108904 


-11 


2.6 


0.8 


6s6d 


'D 2 


115166 


101238 


118678 


116194 


-12 


3.1 


0.9 


6s6d 




116152 


103334 


119000 


116857 


-11 


2.5 


0.6 


6s6d 


3 D 2 


116436 


103555 


119339 


117284 


-11 


2.5 


0.7 


6s6d 


3 D 3 


116831 


103911 


119688 


117758 


-11 


2.5 


0.8 


6p 2 3 


Po 


117408 


108495 


120875 


118450 


-8 


3.0 


0.9 


6p 23 


Pi 


125338 


114961 


129401 


126440 


-8 


3.2 


0.9 


6p 2 3 


Pi 


128817 


117721 


132754 


129839 


-9 


3.1 


0.8 


6s8s 


3 Si 


133568 


120147 


136369 


134187 


-10 


2.1 


0.5 


6s8s 


'So 


134292 


121089 


137132 


134950 


-10 


2.1 


0.5 


6s6p 


3 Po 


49451 


41719 


52320 


50288 


-16 


5.8 


1.7 


6s6p 


3 Pi 


52394 


44743 


55114 


53060 


-15 


5.2 


1.3 


6s6p 


3 P2 


61728 


61728 


65044 


62669 


-14 


5.4 


1.5 


6s6p 


'Pi 


75663 


75663 


76866 


76145 


-7 


1.6 


0.6 


6s7p 


3 Po 


119361 


119361 


122299 


120155 


-11 


2.5 


0.7 


6s7p 


3 Pl 


119576 


119576 


122602 


120472 


-11 


2.5 


0.8 


6s 7p 


3 P2 


122209 


122029 


124873 


122675 


-11 


2.3 


0.5 


6s7p 


'Pi 


122379 


122379 


126014 


124019 


-9 


3.0 


1.3 


6s5/ 


3 F 2 


136216 


136216 


138873 


136600 


-10 


2.0 


0.3 


6s5/ 


3 F 3 


136115 


136115 


138868 


136577 


-10 


2.0 


0.4 


6s5/ 


3 F 4 


136230 


136230 


138870 


136595 


-10 


1.9 


0.3 


6s5/ 


'F S 


136263 


136263 


138997 


136756 


-10 


2.0 


0.4 



TABLE II: Comparison between experimental [13 . [3] and 
Cl+all-order transition energies in cm" 1 . The relative differ- 
ences are given in the last column in percent. 





Transition 


Expt. 


Cl+all-order 


Dif. (%) 


6s6p 


3 P - 6s7s 3 Si 


55778 


55739 


0.07% 


6s6p 


3 Po - 6s6d 3 Di 


66701 


66569 


0.20% 


6s6p 


3 P - 6p 2 3 Pi 


75887 


76152 


-0.35% 


6s6p 


3 Po - 6s8s 3 Si 


84117 


83899 


0.26% 



valence electrons. The ionic core polarizability is evalu- 
ated in the random-phase approximation (RPA), an ap- 
proach that is expected to provide core values accurate 
to better than 5% |4j. We approximate the vc term by 
adding vc contributions from the individual electrons, i.e. 
a„ c (6s 2 ) = 2a„ c (6s), and a vc (6s6p) = a vc (6s) + a vc (6p). 
For consistency, this term is also calculated in RPA. We 
note that a vc contributions are small, but their contri- 
bution to the Aev( 3 .Po — 1 Sq) polarizability difference is 
significant, 15%, due to severe cancelation of the valence 
polarizabilities of these two states. The valence polariz- 



ability is determined by solving the inhomogeneous equa- 
tion of perturbation theory in the valence space, which 
is approximated as 

(E v - H cS )\V(v, M')) = D eff , g |* («, J> M)> (2) 

for a state v with the total angular momentum J and 
projection M [l5jj. The wave function ^(v, M'), where 
A f = M + q, is composed of parts that have angular 
momenta of J' — J, J±l from which the scalar and tensor 
polarizability of the state \v,J,M) can be determined 
[15j [. The effective dipole operator D e ff includes RPA 
corrections. 

Unless stated otherwise, we use atomic units (a.u.) for 
all matrix elements and polarizabilities throughout this 
paper: the numerical values of the elementary charge, 
e, the reduced Planck constant, h — h/2n, and the 
electron mass, m e , are set equal to 1. The atomic 
unit for polarizability can be converted to SI units via 
a/h [Hz/(V/m) 2 ]=2.48832xl0" 8 a (a.u.). The conver- 
sion coefficient is 4ireoa,Q/h in SI units and the Planck 
constant h is factored out in order to provide direct con- 
version into frequency units; ao is the Bohr radius and eo 
is the electric constant. 

While we do not use the sum-over-state approach in 
the calculation of the polarizabilities, it is useful to es- 
tablish which levels give the dominant contributions. We 
evaluate several leading contributions to polarizabilities 
by combining our values of the El matrix elements and 
energies according to the sum-over-states formula for the 
valence polarizability Q: 



\(v\\D\\n)\ 2 
3(2 J 4- 1) ^ E n -E v 



E 



(3) 



where J is the total angular momentum of state v, D is 
the electric dipole operator, and Ei is the energy of the 
state i. 

The breakdown of the contributions to the 6s 2 1 S'o 
and 6s6p 3 Pq polarizabilities cvo of Tl + in a.u. is given 
in Table IIIII Absolute values of the corresponding re- 
duced electric-dipole matrix elements are listed in col- 
umn labeled "D" in a^e. To demonstrate the size of the 
correlation corrections, we list valence results obtained 
in the CI, CI+MBPT, and Cl+all-order approximations. 
The contribution of the other terms listed in the row 
"Other" is obtained by subtracting the sum of the con- 
tributions that are calculated separately from the total 
valence polarizability result obtained by the direct solu- 
tion of the Eq. @ . With the exception of the last column 
labeled CI+all B , we use the theoretical energies obtained 
in the respective approximations. To obtain data listed 
in the last column, we combine Cl+all-order El matrix 
elements and experimental energies. The polarizability of 
the ground state changes by 0.5% as expected from the 
accuracy of the 6s 2 1 Sq— 6s6p Pi transition energy listed 
in TablelH The polarizability of the excited 6s6p 3 Pq state 
remains the same to four significant figures. Such re- 
markable agreement is due to excellent, 0.07%, accuracy 
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TABLE III: Contributions to the 6s 2 1 So and 6s6p 3 Po polarizabilities in a.u. The dominant contributions to the valence 
polarizabilities are listed separately with the corresponding El matrix elements given in columns labeled D. The remaining 
valence contribution is given in row Other. The contribution from the core and vc terms are given by a c and a vc , respectively. 
The dominant contributions to ao listed in columns CI+all A and Cl+all 8 are calculated with CI + all-order energies and NIST 
[l3l. [l4| energies, respectively. The differences of the 3 Po and 1 So polarizabilities calculated in different approximations are 
given in the last row. 



State 



Contribution 



CI 



CI+MBPT 



CI+alE 



ci+air 









D 


ao 


D 


Qo 


D 


ao 


QO 


6s 2 'So 


6s 2 1 S - 


6s6p 3 Pi 


0.424 


0.589 


0.658 


1.149 


0.597 


0.984 


0.997 




6s 2 'S - 


6s6p 'Pi 


2.789 


16.131 


2.619 


13.057 


2.646 


13.450 


13.535 




Other 






0.269 




0.143 




0.155 


0.155 




etc 






4.983 




4.983 




4.983 


4.983 










-0.071 




-0.071 




-0.071 


-0.071 




Total 






21.901 




19.261 




19.501 


19.599 


6s6p 3 P 


6s6p 3 P 


- 6s7s 3 Si 


1.044 


3.113 


0.975 


2.499 


0.980 


2.519 


2.517 




6s6p 3 P 


- 6s6d 3 Di 


2.007 


9.563 


1.893 


7.860 


1.897 


7.912 


7.897 




6s6p 3 P 


- 6p 2 3 Pi 


1.616 


5.219 


1.557 


4.603 


1.562 


4.690 


4.706 




Other 






1.782 




1.630 




1.660 


1.660 




Ctc 






4.983 




4.983 




4.983 


4.983 










-0.338 




-0.338 




-0.338 


-0.338 




Total 






24.322 




21.236 




21.426 


21.425 


Aa ( 3 P - 'So) 








2.421 




1.975 




1.925 


1.826 



TABLE IV: Contributions to dynamic corrections r\ for 
6s 2 'So and 6s6p 3 Po states. 



State 


Transition 


Vn 


n 


?7(6 S 2 'So) 


6s 2 'So - 6s6p 3 Pi 
6s 2 'So - 6s6p 'Pi 


363 
251 


0.000099 
0.000015 
0.000114 


V {6s6p 3 P ) 


6s6p 3 P - 6s7s 3 Si 
6s6p 3 P - 6s6d 3 Di 
6s6p 3 P - 6p 2 3 Pi 


268 
320 
364 


0.000031 
0.000068 
0.000031 
0.000130 


Ar?( 3 P - 'So) 






0.000016 



where y n — u) nv /T; a Q is the static dipole polarizabil- 
ity of the state v, and J the total angular momentum 
of the state v. We list the dominant contributions to r\ 
of the clock states calculated using CI+all-order El ma- 
trix elements and experimental energies in Table |IV] The 
sum in the expression for r\ above converges very rapidly 
making all other contributions negligible. The values of 
r\ for the 6s 1 Sq and 6s6p 3 Pq state are almost equal, and 
their difference listed in the last row gives only 0.0016% 
contribution to the BBR shift. 



IV. EVALUATION OF THE UNCERTAINTY 
AND CONCLUSION 



of the CI+all-order 6s6p 3 Pb — 6s7s 3 Si transition energy 
and opposite signs of the difference between CI+all-order 
6s6p 3 Po—6s6d 3 Di and 6s6p 3 Po—6p 2 3 P± transition en- 
ergies and experiment (see Table HI]). We note that while 
the change in the ground state polarizability is only 0.5%, 
the corresponding change in the final polarizability dif- 
ference Aa( 3 P - 1 So) is 5%. 

We have also calculated the dynamic correction r\ 
of both clock states. The total dynamic correction r\ 
in Eq. ([T]) is the difference of individual corrections, 
At7( 3 P - 1 S ) = v( 3 Po) ~ vC~So). The dynamic cor- 
rection r\ of the state v is evaluated using the formula § 



y, (80/63)tt 2 \(v\\D\\n}\< 



n 



a T (2J+ l)y 3 



21tt 2 

5 Vn 



336tt 4 



We use Table IIIII to evaluate the uncertainty to the 
BBR shift due to the core- valence correlation corrections 
by comparing the CI, CI+MBPT, and CI+all-order re- 
sults for Aa ( 3 P - 1 S ) listed in the last row of Table Iffll 
The difference between the CI and CI+MBPT results is 
23%, which is expected owing to poor agreement of CI 
energies with experiment. The difference between the 
CI+MBPT and CI+all-order results is only 3%. As we 
noted above, the use of the experimental energies changes 
CI+all-order value by 5%. 

We studied the effect of the Breit interaction by re- 
peating the CI+all-order calculation with the one-body 
part of the Breit interaction incorporated into the DF 
equations and construction of the basis set on the same 
footing with the Coulomb interaction. We find that the 
Breit interaction affects both 1 S'o and 3 Pq polarizabilities 
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TABLE V: BBR shifts at T = 300if in B+, A1+, In+, and T1+. B+, A1+, and In+ values are taken from Ref. 0. Polarizabilities 
ao and their differences Aao are given in a.u.; clock frequencies vq and the BBR shifts A^bbr are given in Hz. Uncertainties 
in the values of A^bbr/^o are given in column labeled "Uncertainty". 



Ion 


c*o(%) 


Qo( 3 Po) 


A« 


A^BBR (Hz) 


v (Hz) 


|A^BBR./fo| 


Uncertainty 


B+ 


9.624 


7.772 


-1.85(19) 


0.0159(16) 


1.119 x 10 lb 


1.42 x 10 -1 ' 


1 x 10~ 18 


A1+ 


24.048 


24.543 


0.495(50) 


-0.00426(43) 


1.121 x 10 15 


3.8 x 10~ 18 


4 x 10 -19 


In+ 


24.01 


26.02 


2.01(20) 


-0.0173(17) 


1.267 x 10 15 


1.36 x 10~ 17 


1 x 10 -18 


T1+ 


19.60 


21.43 


1.83(18) 


-0.0157(16) 


1.483 x 10 15 


1.06 x 10" 17 


1 x 10 -18 



by approximately the same amount, —0.5%. As a result, 
the correction to the BBR shift due to the Breit interac- 
tion is negligible (0.6%) at the present level of accuracy. 

To evaluate the uncertainty in the a vc contribution to 
the polarizability, we calculate this term in both DF and 
RPA approximations. The difference between these re- 
sults is taken to be the uncertainty. We find that the 
uncertainty of the vc term contributes 2.4% to the un- 
certainty in the BBR shift. The ionic core polarizability 
a c is the same for both states and does not contribute to 
the BBR shift. 

Based on the comparison of the CI, CI+MBPT, and 
Cl+all-order data, estimated of the accuracy of the a vc 
terms, and estimated effect of the Breit interaction, 
we place an upper bound on the uncertainty of our 
Aa ( 3 -Po — 1 So) polarizability difference and the corre- 
sponding BBR shift of T1+ at 10%. 

Our final result for the BBR shift of the 
6s 2 1 So — 6s6p 3 Po transition in T1+ is A^bbr = 
-0.0157(16) Hz at 300 K. The corresponding relative 
BBR shift at 300 K is |A^ B br/^o = 1-1(1) x 1CT 17 . 
Our final results are summarized in Table El where we 
list the clock state polarizabilities, their difference Aao, 
BBR shift at T = 300 K, 1 S - 3 Pq clock frequencies 
2/q, absolute value of the relative BBR shift IAi/bbr/^oI, 



and the uncertainty in the relative BBR shift of Tl + . 
The Tl + values are compared with the results obtained 
for B + , Al + , and In + ions in Ref. @. The results listed 
in Table [V] demonstrate that near cancelation of the 
1 Sq and 3 Po state polarizabilities in divalent B + , Al + , 
In + ions of group IIIB Q continues for much heavier 
Tl + , leading to anomalously small BBR shift for this 
system. This calculation demonstrates that the BBR 
contribution to the fractional frequency uncertainty of 
the T1+ frequency standard at 300 K is 1 x 10" 18 . We 
find that Tl + has the smallest fractional BBR shift 
among all present or proposed frequency standards with 
the exception of Al + . 
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